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FAR INFRARED EMISSION FROM ELLIPTICAL GALAXIES: NGC 4649, NGC 4472, AND 

NGC 4636 
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ABSTRACT 

We present ISOPHOT P32 oversampled maps and P37/39 
sparse maps, of three bright elliptical galaxies in the Virgo 
Cluster. The maps reach the limiting sensitivity of the 
ISOPHOT instrument at 60, 100, 170 and 200/im. Two el- 
liptical galaxies show no emission at all far-IR ISOPHOT 
wavelengths at a level of few tens of mJy. The null detec- 
tion provides a test of the evolution of dust in elliptical 
galaxies and its size distribution and composition. As pre- 
vious studies have shown, in many elliptical galaxies both 
IRAS and ISO have detected mid-IR excess 6-15 micron 
emission relative to the stellar continuum indicating emis- 
sion from circumstcllar dust. Under the assumption that 
these dusty outflows from evolving red giant stars and 
planetary nebulae are continuously supplying dust to the 
interstellar medium, we have computed the infrared lumi- 
nosity at the ISOPHOT bands appropriate for NGC4472. 
The null far-IR ISOPHOT observations exceed the far-IR 
flux expected from dust expelled from a normal old stellar 
population. 
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1. INTRODUCTION 

In recent years the traditional view of elliptical galaxies as 
simple systems of non interacting stars, devoid of interstel- 
lar matter, has radically changed. Observations across the 
electromagnetic spectrum have demonstrated that the in- 
terstellar medium (ISM) in elliptical galaxies contains sub- 
stantial amounts of cold gas and dust in addition to hot 
gas, the dominant component. Far-IR emission detected 
by the IRAS satellite and improved optical imaging of el- 
liptical galaxies provide clear evidence for the presence of 
dust. However, estimates of the total amount of dust, as 
well as its origin and spatial distribution, remain uncer- 
tain and controversial. Optical observations indicate dust 
masses that are one order of magnitude less than those in- 
ferred by IRAS. We have recently embarked on a program 
to study the infrared emission from early-type galaxies us- 
ing the large database of observations taken by the ISO 
satellite. Our principal goals are to observe the spatial lo- 
cation and emission spectrum of dust throughout elliptical 



galaxies, and with this information, to determine the ori- 
gin, evolution, and physical properties of dust in massive 
elliptical galaxies. 



2. OBSERVATIONS AND DATA REDUCTION 

ISO provides a vast amount of data with good spatial reso- 
lution in the mid and far-IR for a large number of elliptical 
galaxies. A good fraction of these observations have been 
taken using ISOPHOT in several observing modes (Astro- 
nomical Observing Template), including the oversampled 
maps P32, the P37/38/39 sparse maps and P22 multi- 
filter photometry data. Here we present data from three 
very bright ellipticals located in Virgo Cluster: NGC4649, 
NGC4472, and NGC4636. 

NGC4636 has been observed in the P37/39 AOT mode at 
60, 100, and 180 /im. Observations were made in each fil- 
ter with one on-source single point staring mode and one 
off-source exposure. The background position was located 
^ 10' north-east from the target, on blank sky position. 
Data reduction and calibration were performed with the 



PIA 9.1 package (Gabriel et al. 1997). The reduction in- 
cluded correction for the non linear response of the detec- 
tors, readout deglitching, and linear fltting of the signal 
ramps. After resetting of all ramp slopes and subtract- 
ing the dark current, the flux densities were extracted 
using the power calibration of the reference lamps. Af- 
ter sky subtraction the on-source signal was corrected for 
the fraction of the point-spread function included within 
the detector's field of view. A color correction was applied 
assuming that the SED of the galaxy could be approxi- 
mated by a blackbody curve with temperature of 30 K. We 
tested the ClOO detector response on a number of galaxies 
observed in the same observing mode as NGC4636. Inte- 
grated fluxes from 8 ellipticals were compared to the IRAS 
values at 60 and 100 //m. Calibrated survey scans from the 
IRAS satellite were extracted using the SCANPI software 
from IPAC. A very good linear correlation between the 
ISO and IRAS fluxes is seen at both wavelengths. 
NGC4649 and NGC4472 have been observed in P32 mode 
in three broadband filters at 60, 90 and 180 /im using 
the ClOO and C200 detectors. The maps were obtained 
by scanning the spacecraft Y and Z axes in a grid of 4x4 
(NGC4649) and 5x4 (NGC4472) points. The focal plane 
chopper was stepped at intervals of one-third of the detec- 
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Figure 1. NGC4649 P32 maps at 60, 90, and 180 fj,m respectively with optical contours superimposed. The galaxy is 
undetected at all wavelengths observed with ISOPHOT. The bright source at the edge of map is its spiral companion 
NGC4647 and it is prominent at all three wavelengths. The angular resolution of the oversampled map nicely resolves 
NGC4647 and shows that there is no contribution from NGC4649 in this aperture centered on the elliptical. At 180 
fim the signal from the spiral galaxy extends almost to the location of NGC4649 nucleus, but an analysis of the light 
profile shows that the emission is clearly coming from the extended far- IR image of NGC4647. 



tor pixel, at each raster position, providing a sky sampling 
in the Y direction of 15" and 30" for the ClOO and C200 
detectors. To reduce the P32 data we used a dedicated 
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Figure 2. Optical image ( contours ) superimposed to 
ISOPHOT gOtim map for NGG4472 . Similar maps at 60 
and ISO^m show no detectable emission from the bright 
elliptical. 



software package developed at MPI Kernphysik in Hei- 
delberg and supported by the VILSPA ISO Data Center 
( fTuffs fc Gabriel 2002| ). The new routines allow a proper 
correction for transients in PHT32 measurements. Both 
NGC4472 and NGC4649 show no detectable emission at 
60, 90, and 180 /im. The maps reach the limiting sensitiv- 
ity of the ISOPHOT instrument at a level of few tens of 
mJy and the 3cr upper limits are presented in Table 1. Fig- 
ure 1 shows the 60, 90, and 180/im maps for NGC4649. 
The bright spiral companion NGC4647, recorded at the 
edge of the maps, is well detected at all three wavelengths 
and is responsible for the high flux detection, erroneously 
attributed to NGC4649, reported in literature from IRAS 
observations. In fact the angular resolution and beam size 
of the IRAS instruments were not able to disentangle the 
contribution of the two galaxies. Figure 2 show the 90 /im 
oversampled map for NGC4472. 



Table 1. ISOPHOT flux densities. 









Flux (mJy) 






NGC4649 


NGC4472 


NGC4636 




< 137" 


< 147^= 


187± 57 




< 85" 


< 99^= 


491± 64 


Fiao^m 


< 110'= 


< 87" 


790± 71 



" No detection, 3 a upper limit flux 



3 



3. A MODEL FOR NGC4472 



Mid-IR (6 - 15/im) ISO observations of NGC4472 detect 
strong dust emission that can be understood as emission 
from circumstellar dust from mass-losing AGB stars at 
or near the globally expected rate for this galaxy (Athej; 
et al. 2002). Since this dust should emit in the far-IR 



when it becomes interstellar, it may be significant that 
NGC4472 is not detected at A > 25/xm either by ISO 
or IRAS (Knapp et al. 1989). We now briefly describe 



an estimate of the far infrared emission from NGC4472, 
assuming that the grains are dispersed in the hot inter- 
stellar gas before sputtering commences. First we calcu- 
late the temperature Td{r,a) of dust grains of radius a 
(in /Ltm) at galactic radius r and then determine the grain 
emissivity and luminosity for a distribution of initial grain 
sizes. Since the 9.7/xm silicate feature is seen in emission 
in NGC4472, we assume that al l grains have proper ties 
similar to astronomical silicates (Laor &: Draine 1993). 



The grain temperature is determined by both absorption 
of starlight and by electron-grain collisions: 
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47r J*(r, X)Qabs{a, \)'Ka^d\ + 47ra^(l/4)ne(we-Ee)r(a) 



= Ana'^asBTdir, a)^{Qabs){Td,a) 



(1) 



although we consider each type of heating separately. The 
mean intensity of starlight at wavelength A (in fim) is re- 
lated to the B-band intensity J*(r, A) = JxB{r)4){\) and 
the SED cf){\) taken from fTsai fc Mathews (1995)| . JxB{r) 
is found by integrating over a de Vaucouleurs stellar profile 
with a "nuker" core of slope p* cx r~°-^^ within 0.031i?e 
( |Faber et al. 19971 ). At a distance d = 17 Mpc NGC4472 
has an effective radius Re = 8.57 kpc and total stellar 
mass M^.t = 7.26 x 10" M© assuming AU^t/LB = 9.2 



Figure 3. Temperature of dust grains due to starlight heat- 
ing only (dashed lines), and due to electron- grain colli- 
sions only (heavy solid line for grain size a = 1 fj,m; thin 
solid line for grain size a — 0.1 ^m). For grains with size 
a < 0.3 /im the collisional heating is larger than the ra- 
diative one. 



where s = 3.5 and a is in //m (Mathis et al. 1977), then 
the steady state grain size distribution in the hot gas is 



N{r, a) 



da So{r) 
H (s - 1)' 
The coefficient 



a < Qr, 



47rpgl0-i2 



( |van der Marel 1991]) . The second term on the left repre- 
sents heating by e^-grain collisions including a correction 
r(a) for s mall grains which do not completely stop the 
electrons ( Dwek 1986 ). For astronomical silicate grains 
of radius less than ~ 1/im we can use the approxima- 
tions Qabs ~ a'0(A) and (Qabs) ~ 1-35 x 10~^T|a. For 
starhght heating the grain temperature Td{r) — 44JjJ^^ 
K is independent of a and varies from 26 K at 1 kpc 
to 16 K at 10 kpc. For e~-grain heating the grain tem- 
perature Td = 0.662(neT3/2r/a)i/'^ K varies as a-^/^ for 'dt " "/'^p3.2xl0 
a > 0.04/im where r « 1 but is independent of grain 
radius for a < 0.04/im where t « 33a. Figure 3 shows 
the calculated dust temperature profile for the starlight 
heatin g and the o~ grain heating. For NGC 4472 the col 



(2) 



(3) 



depends on the rate that the old stellar population ejects 
mass, a, = 4.7 x 10"^" s~^, the initial dust to gas mass 
ratio scaled to the stellar metal abundance in NGC4472, 
6 = (l/150)zm(r) and the density of silicate grains, pg = 
3.3 gm cm^"^. The grain sputtering rate 
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da 



2 X 10" K 
f 



2.5 



1 



/ims ^(4) 



is a fit to the rate of Draine & Salpeter (1979) with an 
additional correction / = f{a,T) > 1 to account for the 
enhanced (isotropic) sputtering from small grains (Jurac 



lisional heating dominatoe over the radiative heating for 
^ 0.3 



gram s pcs a 



' .6 fim. 



et al. 1998). The hot gas temperature T and proton den- 
sity Uj, for NGC4472 are taken from Brighenti & Mathews 



If grains are expelled from the stars with an miiial MRN 
size distribution 



S{a,r) = So{r)a' 



(1997) 



The infrared emissivity from grains at galactic radius r is 



daN{r, a)ATra'^10^^Qabs{a, A) 
XTTB{r, a, A) 



(5) 
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where B{Td{r), a, A) is the Planck function. 
The total galactic luminosity at wavelength A is 

Lx = J 47rj(r, A)(4/3)^d(r3) 

and the observed flux at A^ is 



F,. = 10' 



mjy. 



(6) 



(7) 



A-ird'^ c 

Figure 4 compares the computed far-IR spectra for the 
cases of stellar radiation heating and coUisional heating 
with the three ISO flux upper limits. The total IR emission 
(not shown in the figure) is the sum of the contribution 
of the two heating sources for each value of amax- Grains 
with an initial MRN distribution having Umax < 1/im are 
consistent with the null observations of NGC4472 in the 
far-IR. Grains with radii less than about 0.003/im expe- 
rience stochastic temperature excursions, but this com- 
plication should not greatly change our estimated flux in 
Figure 4. 
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starlight heating 
e-collisions heating 
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becomes interstellar and is exposed to the mean galactic 
starlight and to bombardment by thermal electrons in the 
hot ^ 1 keV gas. Grains receive somewhat more energy 
from electrons than from stellar photons. The calculated 
far-IR spectrum peaks at A 180 /im, a spectral region 
not accessible to IRAS, but well covered by ISO. How- 
ever, the computed fluxes are below the observed upper 
limits if a < 0.5 /im. Thus the ISO upper limits do not 
strongly constrain the properties of dust in NGC4472. We 
are currently developing a mo del for the far-IR emission 
for NGC4636 (|Temi ct al. 2002|) . This galaxy shows higher 
far-IR fluxes (Table 1) and should allow us to put strong 
constraints on the origin and the evolution of the intcstel- 
lar dust. 

We thank NASA and NSF for grants that funded this 
research. 
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Figure 4- Computed Far-IR flux in mJy from NGC4472 
calculated for Omax — IfJ-m (heavy lines) and Omax — 
O.lfim (light lines). Solid lines show the grain emission 
based only on electron- grain heating and dashed lines show 
the emission from grains that are heated only by stellar ra- 
diation. The three ISO upper limits are shown with open 
triangles. 



4. RESULTS 



We have assumed that the c ircumstellar dust observed in 
the mid-IR from NGC4472 ( |Athey et al. 2002| ) ultimately 



